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Fermi surfaces and orbital polarization in superconducting CeO 0 .5F 0 .5BiS 2 
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We have investigated the electronic structure of BiS 2 -based CeOo. 5 Fo. 5 BiS 2 superconductor us¬ 
ing polarization-dependent angle-resolved photoemission spectroscopy (ARPES), and succeeded in 
elucidating the orbital characters on the Fermi surfaces. In the rectangular Fermi pockets around X 
point, the straight portion parallel to the k y direction is dominated by Bi 6p x character. The orbital 
polarization indicates the underlying quasi-one-dimensional electronic structure of the BiS 2 system. 
Moreover, distortions on tetragonally aligned Bi could give rise to the band Jahn-Teller effect. 

PACS numbers: 74.25.Jb, 74.70.Xa, 78.70.Dm, 71.28.+d 


Since the discovery of BiS 2 -based superconductors by 
Mizuguchi et al 1 in 2012, various BiS 2 -based systems 
have been discovered including RE(0,F)BiS2 system 
(RE=rare earth element). In a typical RE(0,F)BiS2 sys¬ 
tem, electronically active BiS 2 layers (or BiS planes) are 
sandwiched by RE(0,F) layers, and the F substitution is 
believed to be electron-doping to the electronically active 
BiS 2 layer. In the various RE(0,F)BiS2 systems, numer¬ 
ous theoretical and experimental researches have been 
performed in order to understand their electronic states 
including the pairing mechanism 2 19 . Among them, the 
Fermi surfaces and the band dispersions of RE(0,F)BiS2 
have been investigated by means of angle-resolved pho¬ 
toemission spectroscopy (ARPES) 20 22 . The observed 
non-dispersive k x -k z Fermi surfaces 20 due to the two- 
dimensional-like crystal structure are consistent with the 
decoupled BiS plane with block layer by the F-doping 
reported by local structural studies 23 . This indicates the 
absence of Bi 6p z electrons near Fermi level (Ep). The 
observed band structure and Fermi surface geometry are 
also basically consistent with the band structure calcula¬ 
tions, except for that the area of the Fermi surface is not 
fully reproduced. For instance, in the optimal F-doped 
region (x =0.18), the total area of the Fermi surfaces is 
consistent with the amount of F substitution in a similar 
selenide superconductor La(0,F)BiSe2 system 21 . How¬ 
ever, in the higher doped region {x ~ 0.5), the total 
area of the observed Fermi surfaces was much smaller 
than that expected from the nominal amount of F sub¬ 
stitution in Nd(0,F)BiS2 20,22 . The discrepancy in the 
Nd(0,F)BiS2 system suggests that a part of electrons in¬ 
troduced by F substitution would be localized and do not 


contribute to the Fermi surfaces. 

According to the theoretical studies, in a typical 
RE(0,F)BiS2 system, all the Bi 6p orbitals have almost 
the same band center (E = Ep + 2.9 eV), and the 
Fermi surfaces are constructed from the p x and p y or¬ 
bitals due to their much wider band widths than that 
of p^ 24,25 . Near X point where Fermi surfaces exist, 
the p x and p y orbitals may not hybridize and form the 
straight Fermi surfaces along the k y and k x directions, 
respectively. In this situation, the BiS 2 system can be 
viewed as a quasi-one-dimensional electron system 24 de¬ 
spite its two-dimensional-like layered crystal structure. 
The quasi-one-dimensional character of the electronic 
states tends to enhance the Peierls instability due to 
electron-lattice interaction and/or or the charge fluctua¬ 
tion due to electron-electron interaction and may cause 
the deviation from the band structures obtained by the 
local density approximation (LDA) 26 . 

In this paper, in order to address these issues men¬ 
tioned above (smaller Fermi surface, polaron, quasi- 
one-dimensionality), we have performed polarization- 
dependent high resolution ARPES measurement and 
band-structure calculation on CeOo.sFo.sB^. 

High-quality single crystal of CeOo.sFo.sB^ has been 
prepared by CsCl flux method, and its details are 
given elsewhere 27 . The ARPES measurements with lin¬ 
early/circularly polarized photons were performed at the 
undulator beamlines BL-1/BL-9A of the Hiroshima Syn¬ 
chrotron Radiation Center, Hiroshima University. Both 
endstations are equipped with VG SCIENTA R4000 an¬ 
alyzer. The photon energy was set to be 30 eV for all the 
polarizations. We cleaved the single crystalline sample 
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FIG. 1. (Color online) Schematic diagrams of measurement 
geometries with linearly polarized photons that give us p- 
polarization and s-polarization. The corresponding polariza¬ 
tion vectors are denoted by i p and £ s .The p-orbitals with even 
and odd symmetries with respect to the mirror plane are also 
shown in the right. 


in situ in ultrahigh vacuum (< 1 x 10 -10 Torr) in order 
to obtain clean (001) surface. The total energy resolu¬ 
tion including both monochromator and electron energy 
analyzer was measured to be 21 meV. All the measure¬ 
ments were taken at 50 K. The experimental p- and s- 
polarization setups are shown in Figs. 1. In photoemis¬ 
sion process, the transition matrix element can be 
described as 


|M*M|<^|e.r |^>|| 2 ( 1 ) 

where (0^) is the initial- (final-) state wave function, 
and i is the polarization vector of the incident photons. 
The final state <pf is always even since it is in the mirror 
plane and can be approximately described by a plane 
wave. In the p- ( s -) polarization setup, the term e • r is 
even (odd), and the corresponding initial state should be 
even (odd) with respect to the mirror plane in order to 
give non-vanishing Therefore, the p x and p z orbitals 

are observable for the p-polarization setup, while the p y 
orbital for the 5-polarization setup, as depicted in Fig. 
1. Unlike the linear polarization, the dipole selection 
rule does not hold for circularly polarized photons giving 
non-zero for all the orbitals. 

Band structures were calculated using the SIESTA 
method 28,29 , implementing the generalized gradient ap¬ 
proximation in the shape of the Perdew, Burke, and 
Ernzerhof functional 30 . It uses norm-conserving pseu¬ 
dopotentials to replace the core electrons, while the va¬ 
lence electrons are described using atomic-like orbitals 
as basis states at the double zeta polarized level. We 
performed GGA±U calculations 31 because of the known 
strong correlation of the 4/ electrons in cerium. As no 
atomic coordinates have been reported yet, we obtained 
a geometry from density functional theory, by relaxing 
the structure (both atoms and cell), starting from the 
structure known for LaOo.sFo.sB^. The band structure 
calculations are further detailed in Ref. 25. 

The Fermi surface map of CeOo.sFo.sB^ was taken 
by changing tilt angle with circularly polarized light as 
shown in Fig. 2(a) whereas that with the p-polarized 
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FIG. 2. (Color online) (a) Fermi surface maps of 
CeOo. 5 Fo. 5 BiS 2 taken with circularly polarized light and (b) 
with p-polarization geometry, integrated within ±40 meV. (c) 
Cut denoted by ‘cutl’ in (a), (d) MDCs, (e) fit on MDCs from 
E — —0.12 eV up to Ff for the extraction of peak positions, 
(f) Cut denoted by ‘cut2’ in (a), (g) MDCs, (h) fit on MDCs 
from E — —0.16 eV up to for the extraction of peak po¬ 
sitions. Circles in (c,d,f,g) are the peak positions of MDCs, 
dots (solid lines, bars) in (e,h) are the experimental MDCs 
(fit, peak positions) near Ef . Gaussian functions giving the 
best fit to the bottommost MDCs are also displayed in (e,h). 
All the MDCs were obtained with 20 meV energy step. 


light was taken by changing the azimuthal angle as shown 
in Fig. 2(b). Intensities of the Fermi surface maps were 
integrated within ±40 meV above and below Ef . In 
the rectangular Fermi pockets with p-polarized light [Fig. 
2(b)], the straight portion parallel to the k y direction is 
much more enhanced than that parallel to the k x direc¬ 
tion. This observation may indicate that the straight 
portion parallel to the k y direction is derived from Bi 
6p x and that the other straight portion parallel to the k x 
direction has more Bi 6p y character. However, since the 
Fermi surface was measured by changing the azimuthal 
angle, the orbital selectivity is not exact for the straight 
portion parallel to the k x direction. Figures 2(c) - 2(h) 
show detailed analyses on the electronic states forming 
the Fermi surfaces around X point. Figures 2(c) and 
2(f) show the band dispersion of cutl and cut2 denoted 
in Fig. 2(a), and corresponding momentum distribution 
curves (MDCs) are displayed in Figs. 2(d), 2(e) and Figs. 
2(g), 2(h), respectively. In order to estimate the peak po¬ 
sitions, multi-peak fit using two or four Gaussian func¬ 
tions has been done, and the results are shown in Figs. 
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2(e) and (h) (in the present system with strong atomic 
disorder 23 , the MDC width may have rather Gaussian 
character). The Gaussian functions giving the best fit 
to the bottommost MDCs are displayed in the figures. 
The peak positions identified by this fit are denoted by 
circles in Figs. 2(c), 2(d), 2(f), and 2(g) and by vertical 
bars in Figs. 2(e) and 2(h). The results clearly indicate 
that there exist two electron-like Fermi surfaces around X 
point. Full-width-half-maximum was fixed for all the fits. 
The Fermi wave numbers are estimated as kpi ~ O.lA -1 , 
kp 2 ~ 0.12A- 1 , and kp 3 00 0.22A -1 . These labels are in¬ 
dicated in the Figs. 2(e) and 2(h). Identifying the Fermi 
wave numbers enables us to estimate the approximate 
amount of electron in the system; the inner and outer 
electron pockets respectively enclose 3.9% and 7.2% of 
the Brillouin zone. Counting the Luttinger volume of 
the Fermi surfaces, we estimate the amount of electron 
to be 0.22 electron per Bi in CeOo.sFo.sB^, that is much 
smaller than the nominal value. Moreover, the topology 
of the Fermi surface at x ^0.5 predicted by LDA calcu¬ 
lation is different from our results, and likely to be closer 
to the Fermi surface at x ^0.25. 

According to the scanning tunneling microscopy study, 
Bi defects have been found on the BiS plane, whose 
amount was approximately 3% of Bi site 32 . Since one 
Bi 3+ defect introduces three holes to the BiS 2 layer, 
the defect of 3% can eliminate 0.09 electrons per Bi. 
However, this is not enough to explain the reduction to 
x = 0.22 from x = 0.5. One possibility is that the F ions 
remain trapped in interstitial sites and the amount of the 
doped electrons is reduced from the nominal F content. 
The electron probe microanalysis has found that the ac¬ 
tual F concentration can deviate largely from the nominal 
value in this system 33 . Therefore, further investigations 
are needed for the precise F characterization. 

Figure 3(a) shows the schematic Fermi surface, and 
we measured ARPES data along M-T-X line as depicted 
by the bold line. Raw data and EDC curvatures of the 
valence band (here, EDC stands for energy distribution 
curve) taken with circularly polarized photons are shown 
in Figs. 3(b) and 3(c), respectively. Figures 3(d) - 3(g) 
and 3(h) - 3(k) show the band dispersions along M-T-X 
taken with p- and s-polarization setup, respectively. Raw 
data are displayed in Figs. 3(d) and 3(h), corresponding 
EDC curvatures are shown in Figs. 3(e) and 3(i), EDCs 
are in Figs. 3(f) and 3(j), and the detailed electron-like 
band is shown in Figs. 3(g) and 3(k). The ab initio 
calculation is overlaid on each curvature plots. One can 
see a flat band located about E = —1.3 eV (between Bi 
6p conduction band and valence bands) due to localized 
Ce 4/ electrons for both polarizations in Figs. 3(d) and 
3(h) as reported by Ce 4d-4/ resonance experiment on 
Ce(0,F)BiS2 34 . Comparing the band dispersion taken 
with s- and p-polarizations, spectral features are clearly 
observed with p-polarization whereas they are hardly 
seen with s-polarization as shown in Figs. 3(d) and 3(h). 
The EDCs integrated within k x = ±0.4 A -1 along T- 
X direction are also displayed in Fig. 3(1), and one can 
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FIG. 3. (Color online) The polarization dependence of the 
band dispersions of CeOo.sFo.sB^. (a) Schematic Fermi sur¬ 
face and the cutting line along M-F-X. (b) Raw data and (c) 
EDC curvature along T-X with circularly polarized light, (d) 
Raw data, (e) EDC curvature, (f) EDCs along M-T-X, and 
(g) detailed dispersion around X along T-X with p-polarized 
light in the k x direction, (h) Raw data, (i) EDC curvature, 
(j) EDCs along MT-X, and (k) detailed dispersion around X 
along T-X with s-polarized light in the k x direction. The ab 
initio calculation for the CeOo. 5 Fo. 5 BiS 2 system is overlaid 
on each curvature as well. (1) Comparison of EDCs of the 
electron-like band around X with p-, s-, circular, and the sum 
of s and p polarization integrated in the range of k x = ±0.4 
A -1 along T-X direction. 


see that the spectral weights of p- and s-polarizations 
are drastically changed. Considering that the cuts along 
T-X were taken in the direction of k x , the results sug¬ 
gest us that the electron-like band near X point in the 
k x direction (z.e., the Fermi surface is parallel to the k y 
direction) is dominantly derived from Bi 6p x orbital. In 
other words, the system has quasi-one-dimensional prop¬ 
erty in at least T-X cut. 

According to the minimal model calculation by Usui, 
Suzuki, and Kuroki 24 , the band along T-X consists of 
well-hybridized Bi px and py (namely, close to pure Bi 
p x or Bi p x ) whereas that along T-M is derived from pure 
Bi px or py as in Figs. 4(a) and 4(b). Here, the X- and 
T-axes can be obtained by rotating the x- and p-axes 
by 45 degrees with respect to z- axis. Fig. 4(c) shows 
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FIG. 4. (Color online) (a) Schematic diagrams of the or¬ 
bital polarization along T-M predicted by Usui, Suzuki, and 
Kuroki 24 and (b) observed orbital polarization along T-X. 
They can be transformed back and forth each other by their 
linear combination, (c) A schematic diagram of the orbital 
distribution proposed by Usui, Suzuki, and Kuroki 24 and the 
observed Fermi surfaces, (d) Schematic diagrams of the Fermi 
surface and corresponding band dispersion of Bi 6p x /6p y . 
Nesting vectors are also shown here, (e) Schematic diagrams 
of orthorhombically distorted Bi sites with a' < a, Fermi sur¬ 
face, and band Jahn-Teller effect on Bi 6p x /6p y . 


the proposed orbital distribution and schematic Fermi 
surfaces. The basis can be transformed by their linear 
combinations as 

I Px) = /={[ Px ) + \ Py )}, | py) = T{| p x) _ lpy)} , (2) 

At the doping level of x = 0.5, the band along T-M 
crosses Ep and forms almost strait Fermi surfaces with 
Px (or py) character. Usui, Suzuki, and Kuroki pro¬ 
posed the Fermi surface nesting between the px Fermi 
surfaces or between the py Fermi surfaces can enhance 
the spin fluctuations for the superconductivity 24 . On the 
other hand, at the doping level of x = 0.25, the band 
along T-M does not cross Ep and only the band along 
the T-X forms Fermi pockets (disconnected Fermi sur¬ 
faces) around X point with almost pure p x or p y charac¬ 
ter. In this sense, our results are partly consistent with 
the previous calculation for x=0.25. Considering the fact 
that the observed Fermi surfaces of the superconducting 
samples are similar to those predicted for x=0.25, the 
mechanism of the superconductivity should be discussed 
based on the disconnected Fermi surfaces. 


Here, we have succeeded in the detailed orbital charac¬ 
terization; in the rectangular Fermi surface around ( 7 r, 0), 
the straight portion parallel to the k y -axis is dominated 
by Bi 6p x as schematically shown in Fig. 4(d). On the 
other hand, the straight portion parallel to the k x -axis 
has more Bi 6p y character as suggested by the Fermi sur¬ 
face map in Fig. 2(b) although the orbital selectivity is 
not exact in this region. The rectangular Fermi pockets 
around X point keeps a peculiar quasi-one-dimensionality 
due to the orbital polarization. This situation is partly 
similar to the quasi-one-dimensional Fermi surfaces with 
the orbital polarization for x = 0.5 24 . However, as in¬ 
dicated by the arrows in Fig. 4(d), possible nesting 
channels between the straight portions of the rectangu¬ 
lar Fermi pockets may have different orbital character, 
indicating spin and orbital fluctuations. In addition, the 
Fermi surface nesting may enhance the electron-lattice 
interaction. The electron-electron interaction (spin and 
orbital fluctuations) and the electron-lattice interaction 
would be one of the reasons that the band dispersions ob¬ 
tained by ARPES showed poor agreements with the ab 
initio calculation. Interestingly, the orbitally-polarized 
rectangular Fermi surfaces can give the system a pos¬ 
sibility of band Jahn-Teller effect. If the Bi lattice is 
distorted from tetragonal to orthorhombic, for example, 
the degeneracy of Bi p x /p y will be split as schematically 
depicted in Figs. 4(d) and 4(e) (where a' < a), resulting 
in the distorted Fermi surface with C 2 symmetry. As a 
result of the band Jahn-Teller effect, the system has a 
better nesting condition than the tetragonal symmetry, 
which further enhances the one-dimensionality. The one 
dimensional Fermi surfaces may provide Peierls instabil¬ 
ity to the system. The combination of band Jahn-Teller 
distortion and Peierls instability would be responsible for 
the monoclinic distortion under high pressure. 

In summary, we have studied the electronic structure 
of CeOo. 5 Fo. 5 BiS 2 by means of polarization dependent 
ARPES. From the circularly polarized photon, the doped 
electron has been estimated using Fermi surface as 0.22 
electron per Bi in CeOo.sFo.sB^, that is much smaller 
than the nominal value of 0.5. From the linearly polar¬ 
ized photons, we have succeeded in the observation of or- 
bitally polarized Fermi surface. It indicates that the sys¬ 
tem has quasi-one-dimensional property, that could en¬ 
hance the electronic correlation and/or electron-phonon 
coupling. This might be why the band dispersion ob¬ 
tained by ARPES has poor agreements with the ab initio 
calculation. 
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